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Abstract. Accurate calculations for the ground state of the molecular ions He2~'' 
and HeH^+ placed in a strong magnetic field B > 10^ a.u. (~ 2.35 x lO^^G) using 
the Lagrange-mesh method are presented. The Born-Oppenheimer approximation of 
zero order (infinitely massive centers) and the parallel configuration (molecular axis 
parallel to the magnetic field) are considered. Total energies are found with 9-10 
s.d. The obtained results show that the molecular ions Hcj^ and HeH^"'' exist at 
B > 100 a.u. and B > 1000 a.u., respectively, as predicted in [I] while a saddle point 
in the potential curve appears for the first time at i? ~ 80 a.u. and B ~ 740 a.u., 
respectively. 
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1. Introduction 

The neutron star surface is characterized by enormous magnetic fields, typically of 
B ~ 10^^ G (~ 10^ a.u.jj. Even more, for some young neutron stars - magnetars, the 
magnetic field achieves values up to 5 ~ 10^^ G (~ 10^ a. u.). It raises a general question 
about the behavior of matter under such extreme magnetic fields (see for a review e.g. 
[2]). From a practical point of view, because of the intense magnetic fields, the content 
of a neutron star atmosphere became an issue. What kind of chemical species are 
formed there? This question got major importance after 2002 when the Chandra X-ray 
observatory discovered in the spectra of radiation of the isolated neutron star 1E1207.4- 
5209 two absorption features at ~ 0.7keV and ~ 1.4 keV. 

Extensive theoretical studies of traditional and exotic molecular systems in the 
presence of a strong magnetic field have been carried on (for a review see [1] and 
[3]). It was found that for magnetic fields B < 4.414 x 10^'^ G the most bound one- 
electron systems are the exotic molecular ions He2''' and HeH^"*" besides the traditional 
molecular ion H^. With the hypothesis of a neutron star atmosphere composed of 
mixed hydrogen-helium molecules, a model was proposed [4] that could explain the 
two absorption features discovered by the Chandra X-ray observatory. Even more, if 
the magnetic field strength is increased, new exotic finite molecular chains start to be 
present. In particular, for B = 2.1 x 10^ a. u. the potential curve of the exotic system 
Li2''' start to develop a minimum [5]. 

The aim of the present work is to make an accurate study of the two molecular ions 
He^+ and ReR^^ in the presence of a strong magnetic field B < 4.414 x 10^'^ G where 
the relativistic corrections seem to be unimportant. The molecular axis is assumed to 
be aligned parallel to the magnetic field line because this configuration is optimal if 
a magnetic field is sufficiently strong [3l [6]. The Born-Oppenheimer approximation is 
considered where the positively-charged centers are supposed to be infinitely massive. 

Following [5] in the case of one-electron systems the variational method and the 
Lagrange-mesh method are two methods which complement each other. The first one 
has been proven to be very efficient [3] with physically motivated trial functions. On the 
other hand, the Lagrange-mesh method nowadays provides the most accurate results 
for some simple one-electron molecular systems in the presence of a strong magnetic 
field. In practice, the implementation of the Lagrange-mesh method calculations is 
made easier and often feasible with a priori knowledge of the equilibrium distance and 
the corresponding total energy. Such a knowledge is provided by a variational calculus. 
Thus, following the results presented in [1] the Lagrange-mesh method is applied to check 
and further improve the variational results. 

Atomic units m-e = e = h = 1 are used throughout, although the energy is given in 
Rydbergs. 

t 1 a.u. « 2.35 X 10^ G= 2.35 x 10^ T. 
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2. One-electron molecular ion at the Born-Oppenheimer approximation 

Choosing the vector potential corresponding to a constant magnetic field B = (0, 0, B) 
in the symmetric gauge A = ^{—y, x, 0), the Hamiltonian which describes two infinitely 
heavy centers of charges Zi = 2 and Z2 = 1,2 situated along the z axis and one electron 
placed in a uniform constant magnetic field directed along the 2;— axis, B = (0, 0, B) is 
given by 

2Z^ 2Zo 2 Z\ Zo _ ^ B'^ o ^ s 

(see Fig. [1] for the geometrical setting). Here = + and is the 2;-component of 
the electron orbital momentum, which is an integral of motion, L^] = 0. 



• e 




Figure 1. Geometrical setting of tlie ground state Xog for the helium containing 
molecular system of one-electron and two charged centers Z\=2 and Zi = 1,2 placed 
in a magnetic field B in parallel configuration along the z-axis. R is the internuclear 
distance and r\ , r-i are the distances between the electron and the charged centers Z\ 
and Z2, respectively. 



3. The Lagrange- mesh method 

As it was pointed in [5] the Lagrange-mesh method and the variational method naturally 
complement each other. With this in mind, the Lagrange-mesh method [71|8l|9] is applied 
with the aim to confirm and improve the variational results presented in [1]. The most 
adequate coordinate system to describe a two-center molecular system of charges Z\ 
and Z2 is the system of spheroidal coordinates (^, 77, (p) defined as 

(see Figll]) with ^ G (0, 00), r] G ( — 1, 1) and the azimuthal angle G (0, 27r). Since the 
z-projection of the orbital angular momentum commutes with the Hamiltonian ([1]), 
it is replaced by its eigenvalue m (magnetic quantum number). The wave function for 
the ground stated (m = 0) should be nodeless and it can be written as [7] 

^n.=o(r) = ^^^^^o(e,r/). (3) 
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The function ipoi^^yT]) is expanded in the Lagrange basis as 

Mi.v) = Y.Y.^Mi.v) , (4) 

i=l i=l 

where A^^ is the size in the ^-direction and is in the //-direction. An exphcit expression 
of Fij{^,ri) and more details are given in [7]. 

Taking the wave function (jl]) and using the Gaussian quadratures associated to 
each coordinate, the Schrodinger equation for two-charged centers gets the form of mesh 
equations (see |7]) 

y^jTjji'j/ + V{hxi, -qj) SijSi'j' - E Sii'Sjj']ci'j' = . (5) 
i'=i j'=i 

The kinetic energy matrix elements Tiji/j/ are calculated in [7]. The potential V{C,,ri) 
is evaluated at the zeros of the Laguerre L^^^Xi) = (scaled by the dimensionless 
parameter h) and the Legendre polynomials PNniVj) = 0- Finally, the problem of solving 
the Schrodinger equation is reduced to searching eigenvalues of the matrix equation ([5]). 

4. Results 

Results of the Lagrange-mesh calculations for the m = ground state of the molecular 
ions Hea"*" and HeH^"^ placed in strong magnetic fields B > 10^ a.u. in parallel 
configuration are presented in Tables [1] and [21 respectively. Both exotic systems, do 
not exist without the presence of a strong magnetic field. The results, demonstrate 
that the potential energy curve starts to display a well pronounced minimum at finite 
internuclear distance at the magnetic fields Bt^ ~ 10^ and ~ 10^ a.u., for the ions 
He2''' and HeH^"*", respectively. It provides a theoretical indication about the possible 
existence of the bound exotic molecular ions Hef and HeH+. The potential energy 
curve as a function of the internuclear distance R is characterized by the presence of a 
potential barrier. At large internuclear distances the interaction of two charged centers 
becomes repulsive: the total energy curve approaches from above to the total energy of 
the atomic ion He'''. The system is metastable towards the dissociation to He'*' +p{a) if 
the energy of He'*' is lower than the minimum on the potential curve. Both considered 
systems display two characteristic properties of the Coulomb systems in a magnetic field: 
as the strength of the magnetic field increases, they become more bound (the binding 
energy grows) (i) and more compact (the equilibrium distance decreases) (ii). 

Tables [1] and [2] present the equilibrium internuclear distance Rgq, the corresponding 
total E^'^"' and binding Ef, = B — E^"^^ energies, the He+ total energy from which the 
dissociation energy Erfi,, = Er"-^P'' can be obtained, as well as the position Rmax of 
the maximum E^""^ of the barrier and its height i.e. the difference AE = E^"-^ — E^^"' 
for the molecular ions Hej^ and HeH^'*', respectively. 

The lowest vibrational energy E^^'' is calculated by using the harmonic oscillator 
approximation around the equilibrium position. It allows us to define the zero point 
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energy and is presented in Tables [T] and [21 The computer code JADAMILU [lOj was 
used for the lowest-eigenvalue search in the mesh calculations. 

4.I. The molecular ion He2^ 

Table 1. Ground state lag of the molecular ion He2^ in a magnetic field B: 
equilibrium internuclear distance Reg, total i<;™™ and binding Et, = B — E""^'' energies; 
He"*" and binding energies; Rmax of the maximum at the potential curve; height of 
the barrier A£' = iJ™'^^ — and lowest vibrational energy E^^^ . For each value of 

the magnetic field the second line presents the variational results from [I] 



B 


Req 




Et[Ry] 






Rrna. 


AE[Ry] 






10^ a.u. 


0.77926 


83.47593826 


16.52406174 


19.12109305 


10.292162923 


1.0219 


0.0331 


0.026 


Present 


0.780 


83.484 


16.516 


19.11 


10.291 


1.02 


0.033 


0.026 


m 


1012 G 


0.42004 


396.8490192 


28.6828956 


30.93148909 


17.1522314 


0.9044 


1.0250 


0.100 


Present 


0.420 


396.864 


28.668 


30.87 


17.143 


0.90 


1.024 


0.100 


DP 


10^ a.u. 


0.30923 


960.7139916 


39.2860084 


40.5413906 


22.7816779 


0.8140 


2.4633 


0.170 


Present 


0.309 


960.732 


39.268 


40.40 


22.779 


0.82 


2.466 


0.169 


HI 


10"G 


0.19314 


4190.162674 


65.156475 


62.392739 


35.760513 


0.6818 


7.3154 


0.369 


Present 


0.193 


4190.182 


65.137 


61.99 


35.754 


0.70 


7.328 


0.366 


m 


10* a.u. 


0.14955 


9913.740241 


86.259759 


79.097975 


45.80516 


0.6175 


12.2422 


0.544 


Present 


0.150 


9913.767 


86.233 


78.43 


45.797 


0.62 


12.25 


0.561 


m 


4.414x101^0 


0.1256 
0.126 


18677.83406 
18677.857 


105.14466 
105.121 


93.45755 
92.53 


54.54632 
54.502 


0.576 
0.58 


17.1720 
17.19 


0.732 
0.739 


Present 

m 


10'' a.u. 


0.08174 


99828.5786 


171.4214 










1.451 


Present 



The results for the ground state IcXg of the symmetrical system Heg for magnetic 
fields B > 10^ a.u. in parallel configuration are presented in Table [H For each value 
of the magnetic field shown in Table [H the results are compared with those presented 
in Turbiner and Lopez Vieyra [1] in the frame of the variational method. In all cases 
the Lagrange-mesh total energies are essentially more accurate than the corresponding 
variational results. For the binding energies, the relative accuracy is ~ 0.05% at -B = 10^ 
a.u. while it is ~ 0.02% at B = 4.414 x lO^^^ G. 

Increasing the magnetic field value from zero, the potential curve first does not 
present an indication of any bound state. However, for 5 ~ 80 a.u. a saddle 
point appears and eventually for a magnetic field B ~ 10^ a.u. the potential energy 
curve starts to display a well pronounced minimum for R ~ 0.78 a.u. indicating the 
existence of a metastable state of the molecular ion Hcg^, unstable towards decay to 
He2''' He"*" + a. For this value of the magnetic field the energy curve shows that 
the potential barrier is high enough to keep at least one vibrational level, E^''' < AE. 
Hence, the system is stable against vibrations. The height of the barrier increases when 
the magnetic field grows and the system becomes more stable allowing more vibrational 
modes. Both AE and E^^^'' increase monotonously with the magnetic field, their ratio 
AE/E^'^ increases from ~ 1 for 5 = 10^ a.u. up to ~ 23 for B = AAU x 10^^ G. At 
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the same time, the dissociation energy E^iss = — E^'^'^ decreases and eventually, 
for mae; netic fields B ~ 10^^ G., the total energy becomes smaller than the total 
energy of the He"*" atomic ion. Thus, the system becomes stable towards the decay 
He2''' — )• He"*" + a. The qualitative behavior of this molecular system is typical: when 
we increase the magnetic field strength, the molecular ion Hcg'*' becomes more compact 
(the equilibrium distance decreases) and more bound (the binding energy increases), see 
e.g. [3]. 

4-2. The molecular ion HeH^^ 

Table 2. Ground state lag of the molecular ion HeH^^ in a magnetic field B: 
equilibrium internuclear distance R^q, total El"™ and binding Eb = B — E"''™ energies; 
He"*" and binding energies; Rmax of the maximum at the potential curve; height of 
the barrier AE' = iJ™'*^ — and lowest vibrational energy E^'^^. For each value of 

the magnetic field the second line presents the variational results from [I] 



B 




Et'""'[Ry] 
















10'' a.u. 


0.3161 


962.6035199 


37.3964801 


40..5413906 


22.7816779 


0.4215 


0.0581664 


0.126 


Present 


0.316 


962.635 


37.365 


40.40 


22.779 


0.424 


0.061 


0.13 


Q 


IQi^G 


0.1846 


4195.651303 


59.667846 


62.392739 


35.760513 


0.3756 


1.189682 


0.413 


Present 


0.185 


4195.693 


59.626 


61.99 


35.754 


0.375 


1.197 


0.41 


HI 


10* a.u. 


0.14073 


9922.639585 


77.360415 


79.09797 


45.80516 5 


0.3416 


2.65622 


0.604 


Present 


0.142 


9922.697 


77.303 


78.43 


45.797 


0.341 


2.651 


0.65 


m 


4.414xl0i-^G 


0.1190 
0.119 


18690.07669 
18690.121 


92.90203 
92.858 


93.45755 
92.53 


54.54632 
54.502 


0.3185 
0.318 


4.21166 
4.22 


0.878 
0.88 


Present 

C! 



Table [2] presents the results for the la state of the asymmetric two-center system 
HeH2+ for magnetic fields B > 10 a.u. in parallel configuration. For each value of the 
magnetic field shown in Table [21 the lower line corresponds to the variational results 
given in pQ. Systematically, the Lagrange- mesh results give much more accurate results 
for all studied magnetic fields. For the binding energy, the relative improvement of the 
mesh calculations is ~ 0.08%, for B = 10^ a.u. up to - 0.05% for B = 4AU x 10^^ G. 

As a function of the magnetic field a saddle point appears on the potential 
curve when the magnetic field value is i? ~ 740 a.u. Increasing the magnetic field 
until B ~ 10^ a.u., a well pronounced minimum on the potential curve emerges at 
R ~ 0.32 a.u., indicating the formation of a metastable state, unstable for the decay 
to HeH^"*" — > He'*' + p. The potential barrier is small and it does not allow to keep a 
single vibrational level E^^^ > AE. Finally, when the magnetic field strength becomes 
B ~ 10^'^ G, the potential energy well is sufficiently deep to keep more than one 
vibrational state. Also, the height of the barrier increases when the magnetic field 
is increased and the system becomes more stable against vibrations. The ratio AE / E^^^ 
increases from ~ 3 for B = 10^^ G up to ~ 5 for B = 4.414 x 10^'^ G. The dissociation 
energy E^iss = i?^*" — E^^~^ decreases but the system remains unstable towards the 
decay HeH^'*' He^ + p for all magnetic fields studied. The system is always stable 
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towards the decay HeH^"*" — )• H + a due to the fact that the hydrogen atom in a magnetic 
field is the least bound one-electron system. 

5. Conclusions 

It is presented a study in the non-relativistic frame and the Born-Oppenheimer 
approximation of zero order of the ground state of the molecular ions He^+ and ReR^^ 
in the presence of a strong magnetic field B < 4.414 x 10^'^ G in parallel configuration. 
Although the Lagrange-mesh method is applied, the variational results presented in [1] 
are taken into account as a reference point for the mesh calculations (see [5]). 

The results give clear indications that the exotic molecular ions He^+ and ReR^^ 
begin to exist as metastable states starting at the threshold magnetic fields Bth ~ 10^ 
and ~ 10^ a.u., respectively. As the magnetic field increases the potential wells of both 
systems become sufficiently deep to keep more than one vibrational state. Eventually, for 
B = 10^^ G the ion Re^'^ becomes stable. The HeH^^ ion remains unstable towards decay 
HeH^^ — ^ Re'^+p in the whole domain of magnetic fields considered B < 4.414 x 10^^ G. 
For both systems and all studied magnetic fields the obtained results are essentially more 
accurate then those presented in [1]. 

These molecular ions have a similar behavior: when the magnetic field increases, 
each molecular ion becomes more compact and more bound. For all magnetic fields the 
most bound helium-containing molecular system is Hcg^, i.e., E^'^^ > E^'^^^^ . From 
Tables [1] and [2] it can be seen that the molecular system placed in a magnetic field 
is the less bound two-center one-electron molecular system at B > 10^^ G. 
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